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We study the astro-physical impact of inaccurate and incomplete modeling of the gravitational 
waveforms from compact binary coalescences (CBCs) . We do so by the matched filtering of complete 
inspiral-merger-ringdown (IMR) signals with a bank of inspiral-phase templates modeled after the 
3.5 post-Newtonian TaylorTl approximant. The rationale for the choice of the templates is three- 
fold: (1) The inspiral phase of the Phenomenological signals, which are an example of complete IMR 
signals, is modeled on the same TaylorTl approximant. (2) In the low-mass limit, where the merger 
and ringdown phases last much shorter than the inspiral phase, the errors should tend to vanishingly 
small values and, thus, provide an important check on the numerical aspects of our simulations. (3) 
Since the binary black hole (BBH) signals are not yet known for mass-ratios above ten and since 
signals from CBCs involving neutron stars are affected by uncertainties in the knowledge of their 
equation of state, inspiral templates are still in use in searches for those signals. The results from 
our numerical simulations are compared with analytical calculations of the systematic errors using 
the Fisher matrix on the template parameter space. We find that the loss in signal-to-noise ratio 
(SNR) can be as large as 60% even for binary black holes with component masses mi = 13M© 
and 7U2 = 20Mq. Also, the estimated total-mass for the same pair can be off by as much as 20%. 
Both of these are worse for some higher-mass combinations. Even the estimation of the symmetric 
mass-ratio 77 suffers a nearly 20% error for this example, and can be worse than 50% for the mass 
ranges studied here. These errors significantly dominate their statistical counterparts (at a nominal 
SNR of 10). It may, however, be possible to mitigate the loss in SNR by allowing for templates with 
unphysical values of r\. 
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INTRODUCTION 



The next (second) generation of gravitational-wave 
(GW) detectors and their network are being planned to 
come online around 2014. These detectors include the 
LIGO detectors [l[ in Hanford and Livingston, USA, the 
Virgo detector in Pisa, Italy 0, and the GEO600 detec- 
tor in Ruthe, Germany Q. Their quest will not only be 
to make the first detections but also to begin doing as- 
tronomy with them. That entails measuring the source 
parameters as accurately as possible. Measurement ac- 
curacy is critical for triggering searches of electromag- 
netic counterparts in other observatories, which can pro- 
vide more complete knowledge about these astrophysi- 
cal objects by tapping into the complementary channels 
of information accessible to them. That information in 
turn will enrich our understanding of how, where, and at 
what rate these objects are formed, apart from answer- 
ing other questions and raising new ones. This endeavor 
will be limited, on the one hand, by the inherent statis- 
tical noise in the measurement process and, on the other 
hand, by the accuracy with which the search templates 
can model actual gravitational waveforms. The former 
issue of the statistical measurement error was studied 
in detail in Ref. Q for estimating the parameters of 
coalescing binaries of non-spinning black holes. The lat- 
ter issue is one of systematics and is discussed here for 
searches using post-Newtonian templates. 



'Electronic address: sukanta@ws u.edul 
^Electronic address: shaonghosh@wsu.edu 
^Electronic address: ajith@caltcch.edu 



Owing to both the anticipated rates and the theo- 
retical knowledge of their waveforms, coalescing binary 
black holes (BBHs) arc among the most promising GW 
sources sought by ground-based GW detectors. The in- 
spiral and ringdown phases of these waveforms are accu- 
rately known through the post-Newtonian (PN) approx- 
imation to General Relativity and black hole perturba- 
tion theory, respectively. And developments in numerical 
relativity (NR) in the last five years has made it possi- 
ble to compute accurate gravitational waveforms from 
the hitherto unknown merger stage as well 0413 ■ I n 
spite of this progress, our knowledge of the gravitational 
waveform of these and other compact binary coalescences 
(CBCs) is not complete over regions of the parameter 
space that can be astrophysically significant. First, nu- 
merical merger waveforms from BBH systems arc not 
yet available for mass-ratios greater than ten; also, more 
studies are awaited on how robust the published merger 
waveforms for mass-ratios greater than six are. More 
pertinently, complete inspiral-merger-ringdown (IMR) 
waveforms are not yet available for mass-ratios above 
four. Second, for compact binaries involving at least one 
neutron star (NS) the effects of the variety of possible 
NS equations of state (EOS) and magnetic field on their 
waveforms arc still under study. Some progress has been 
made in directly relating the differences among the NS 
EOS to those in the amplitude and phase evolution of 
the waveforms fl3rl2(il ]. But further research is required 
to establish the robustness of these results. Regardless 
of the outcome, unraveling the NS EOS will require a 
high signal-to- noise ratio (SNR), which may not be very 
probable in the era of Advanced LIGO or "AdvLIGO" 
(see, e.g., [2(| and the references therein). However, that 
by itself does not imply that the detectability of these 
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systems will be hampered. Indeed, often the hope is ex- 
pressed that a sufficiently large bank of filters across the 
parameter space of binary masses (and spins), even if 
somewhat inaccurate and incomplete, may still be able 
to capture these waveforms without a significant loss of 
SNR. Nevertheless, it is critical to establish how big or 
small that loss is for aiding the formulation of detection 
strategies for upcoming searches and the planning of fu- 
ture missions. 

To quantitatively assess the effects of inaccurate mod- 
eling, one needs knowledge of the exact waveforms, which 
we do not possess. To break this impasse, we employ a 
strategy that was explored earlier by Cutler and Vallis- 
neri for LISA [2l[ • Namely, we choose a surrogate for the 
exact waveform that, in spite of its approximate nature, 
is useful in estimating the loss in SNR and parameter 
accuracy resulting from our incomplete knowledge of its 
source. 

Specifically, we choose as our surrogate waveforms the 
Phenomenological inspiral-merger-ringdown waveforms 
of Ref. [22]. These are analytic waveforms that have 
been modeled to have better than 99% fitting factor 
with hybrid waveforms constructed from PN and NR 
waveforms, with mass-ratios from one to four. Here, we 
take these Phenomenological waveforms to be the "ex- 
act" waveforms in the extended mass-ratio range from 1 
to 8, which is not much wider than their proven range 
of validity. To assess the drop in SNR suffered from 
an inaccurate modeling of the CBC waveforms, we com- 
pute the fitting-factor of a bank of TaylorTl 3.5PN tem- 
plates [HI when filtering simulated AdvLIGO data with 
an IMR signal in it. In the process, we are also able 
to estimate the error in the maximum likelihood esti- 
mates of both mass parameters arising from inaccurate 
waveform modeling. Similar studies on systematic errors 
were carried out in Refs. [HI, [HI- Those studies differ 
from ours primarily in that they base their target signals 
on the effective-one-body (EOB) model [H||. Also, the 
range of CBC mass pairs for which we numerically com- 
pute the loss of SNR and parameter-estimation biases is 
much wider than those studied in the past. Moreover, 
we vet our numerical results against analytic Fisher ma- 
trix calculations in the parameter region where the two 
are expected to agree. As in Ref. Q, here too we limit 
our attention to the leading harmonic of non-spinning 
CBC signals. The sensitivity we assume for the second 
generation detector is that of Advanced LIGO 27[. 



II. 



BINARY BLACK HOLE SIGNALS AND 
THEIR PARAMETER ESTIMATION 



the ones studied here, these functions can be treated as 
constants in time. 

The two polarization components of the dominant har- 
monic of BBH signals studied here are sinusoids with 
varying amplitude and frequency, and have phases 7r/2 
radians apart relative to each other. Accordingly, their 
GW signal in a detector can be written as: 



h(t) = CA(t) cos[(p(t) + ipo] 



(2.2) 



where the amplitude coefficient C and the initial phase 
po can be taken as constants. The signal's time-varying 
phase tp(t) and amplitude A(t) are functions of the phys- 
ical parameters of the binary, such as the component 
masses (and the spins). 

The interferometric data in which these signals are 
searched for are noisy. Thus, any test for establishing the 
presence of a signal in that data requires the modeling 
of this noise, n(t). We denote its Fourier transform as 



»(/) 



n(t) e~ 2nift dt . 



(2.3) 



and take it to be zero-mean Gaussian and stationary: 



i(t) = 0, 



(2.4) 



n*(f)n(f) = -S h (f)5(f-f). (2.5) 

Above, the over-bar denotes the ensemble average and 
Sh(f) is the Fourier transform of the auto-covariance of 
the detector noise and is termed as its (one-sided) power 
spectral-density. 

We also assume the noise to be additive, which implies 
that when a signal is present in the data s(t), then 



s(t) = h(t) + n(t) 



(2.6) 



The noise covariance Eq. (|2.5p introduces the following 
inner-product in the function space of signals: 



(a, b) = m df 



a *(f)kf) 
S h (f) 



(2.7) 



where a(/) and b(f) are the Fourier transforms of a(t) 
and b(t), respectively. For the above model of the de- 
tector noise the Neyman-Pearson criterion [2^ yields an 
optimal search statistic, which when maximized over the 
amplitude coefficient C, is the cross-correlation of the 
data with a normalized template, 



P = (h,s) 



(2.8) 



The gravitational- wave strain in an interferometric de- 
tector can be expressed in terms of its two linear polar- 
ization components h + (t) and h x (t) as 



h(t) = F+h+{t) + F x h x (t) 



(2.1) 



where F + and F x are the detector's antenna-pattern 
functions. These functions depend on the two angles 
locating the source in the sky and a third angle specify- 
ing the orientation of the polarization ellipse (see, e.g., 
Ref. [HI ). For transient sources lasting not more than 
several minutes in a detector's frequency band, such as 



where the normalized template is h(f) = h(f)/ J (h, h). 
Here, h(f) obeys the condition 



\h\\ 2 = (h, h) = 1 



(2.9) 



and is said to have a unit-norm. When s is replaced by 
its unit-norm counterpart (i.e., by s/||s||) in Eq. (|2.8j) . 
the resulting inner-product is called the match between 
two unit-norm waveforms. 

In a "blind" search in detector data, where none of 
the binary's parameters are known a priori, the search 
for a GW signal requires maximizing p over a "bank" of 
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templates (see, for e.g., [301) corresponding to different 
values of those physical parameters. Apart from the in- 
trinsic source parameters, the waveform also depends on 
the (unknown) initial phase ipo and the time of arrival 
to- Maximization over the initial phase ifo is effected by 
using two orthogonal templates for each combination of 
the physical parameters [3l| , and the maximization over 
to is attained efficiently with the help of the Fast Fourier 



Transform (FFT) algorithms |32|. 

For modeling the signals, we use the analytical Fourier 
domain IMR waveforms proposed in Reference [Hj]: 



h{f) = A eS {f)^M\ 



(2-10) 



where the effective amplitude and phase are expressed 
as: 
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(2.11) 



fe=0 



In the above expressions, 



is 



(/ - /ring) 2 + 



(2.12) 



is a Lorentzian function that has a width a, and that 
is centered around the frequency / r i ng . The normal- 

/ f . \- 2 / 3 

ization constant, w = I ^ lln « j ; i s chosen so as 

to make A c g(f) continuous across the "transition" fre- 
quency /ring- The parameter / mG rg is the frequency at 
which the power-law changes from f~ 7 ^ 6 to /~ 2 / 3 . The 
effective distance to the binary is denoted by d e g, which 
is related to the luminosity distance di, by d e s = dhjC. 
The phenomenological parameters /merg, /ring, c and / cu t 
are given in terms of the total mass M and symmetric 
mass-ratio r\ of the binary as 



s, which is just the discrete counterpart of s(t). Let the 
time-series corresponding to exact and approximate GW 
signals from CBCs, parameterized by 9^, be denoted by 
{hx(# M )} and {1ia(0 m )}, respectively. In general, these 
two sets will lie on two different submanifolds of this 
vector space. Consider the simple case where the lat- 
ter denotes just the inspiral part, sans the merger and 
ringdown parts, of the former. Clearly, the latter will lie 
on a submanifold of the space of the former time-series 
of the complete waveform. If, additionally, the latter 
waveform were different from the former also owing to 
a parameter-dependent normalization factor, then the 
latter can cease to live on a sub-manifold of the former. 

Next, let an exact signal hx(#f r ) be present in the 
detector data, 



7rM/ mcrg 

TT-W/ring 

irMa 
TMf cut 



The coefficients a j ,bj,< 



aoV + bo i] + co , 

ai rf +bl1] + Cl, 



a2V 



b 2 r} 
b 3 i] 



C2 
C3- 



(2.13) 



j = 0...3 and x fc ,y fc ^, k 



For 



0,2,3,4,6,7 are tabulated in Table I of Ref. [33 
component masses m^2, the total mass is M = mi +7712 
and the symmetric mass-ratio is r\ = m\mijM . For the 
discussion here, it helps to remember that for a mass- 
ratio of 7711/7712 = 1, 4, and 8, one has n = 0.25, 0.16, 
and ~ 0.1, respectively. 



A. Parameter measurement errors 

Following Cutler and Vallisneri [2l|, consider the vec- 
tor space of possible time-series data from a GW detector 



s = h x (C) + n. 



(2.14) 



Then the template belonging to the approximate wave- 
form family that gives the best fit to the above data will, 
in general, have parameters 9^ { that arc different from 
6>f r . If we denote that template as {hA(^ rj ' f )}, the best fit 
parameters are those that minimize the distance between 
the signal and the template, namely, 



d 3 \\s-h A (9)\ 



s-h A (0 b f)> =0. 



(2.15) 

For small differences between the waveforms, the above 
expression can be Taylor expanded about bf = 9l T + A9 l 
to show that 



J 



A9 l = (r- 1 ^))^ (^h A (0bf), n) + (d,h A (9 h! ), [h x (0 tr ) - h A (0 tr )] )} 



(2.16) 
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where TpjObi) = (<9ihA(#bf), dj\i^(6bi ))■ As was noted 
in Ref. [21(, whereas the first term above gives the sta- 
tistical contribution to the error in the parameter esti- 
mation, the second term gives the systematic one. The 
relevant term is the one that dominates the other for 
expected signal strengths in a given detector or network. 

A total of nine parameters characterize the non- 
spinning BBH coalescence signals considered here. They 
are the total mass M , the symmetric mass-ratio 77, the 
sky-position angles (a, <5), the binary's orientation angles 
(ip, 1), the luminosity distance cIl, the initial (or some 
reference) phase ipo, and the time of arrival (or some ref- 
erence time) to- Here we present results for systematic 
errors in M and 77, and the fractional loss of SNR, arising 
from inaccurate waveform modeling. Systematic errors 
in the complementary set of parameters, especially, the 
ones obtainable only with multi-site observations, will 
be presented elsewhere HH . 



B. Numerical simulations 

For numerically computing the matched-filter outputs 
and maximizing them over the template parameters, we 
took the complete IMR signals as the exact or target 
signals, with mi. 2 <E [13, 104]/V/q. Thus, the mass-ratio 
of the target signals ranged from 1 to 8. The template 
bank is chosen to comprise 3.5PN TaylorTf waveforms 
[23l | , with mass parameters overcovering the mass-range 
of the target signals. The choice of the template wave- 
forms is governed by the fact that the inspiral phase of 
the Phenomcnological waveforms is modeled after that 
PN approximant. The templates are modeled with M 
and 77 such that 7711,2 € [5,121]M , but always with 
f] < 0.25, which is the physical upper-bound. For these 
studies, we used the method and the code described in 
Ref. [ij . Only one target signal is present in the data at 
any given time. 

Gravitational-wave searches for CBC signals in a sin- 
gle interferometric detector use unit-norm templates 
whereby the signal amplitude is measured from the value 
of the signal-to- noise ratio. The measured values of other 
parameters, such as a binary's component masses, to and 
0o , are those defining the template that yields the maxi- 
mum match with the injected signal. These are the same 
values that maximize the match between a unit-norm 
template and a unit-norm signal. Given a unit-norm sig- 
nal, the fitting factor (FF) of a template bank is defined 
as the maximum match it yields for that signal. The FF 
is a useful construct because (1 — FF) is the fractional 
loss of SNR of the target signal when filtered with the 
chosen template bank, and depends only on the mass 
parameters of the target signals. Specifically, it is inde- 
pendent of the CBC signal parameters that cannot be 
measured by a single detector, namely, the distance and 
the polarization, inclination, and sky-position angles. 

The fitting factor for the above choice of template 
bank and target signal family is presented in Fig. [T] 
where the maximization over the template parameters 
(to , M, rj) was carried out numerically and that over </>o 
was carried out analytically, as explained above. Figure 
[T]rcveals that for any given mass-ratio, the fitting factor 
first dips with increasing M, before recovering somewhat 
for higher values of target M. This behavior is accom- 



panied by an increasingly negative error in AM/M, as 
shown in Fig. [5J and relatively smaller errors in A77/77, 
especially, for the region around target 77 ~ 0.25, as de- 
picted in Fig. |4j This set of observations is explained 
by the fact that, for any target waveform, the templates 
that give the best fit tend to have a smaller M, which 
tends to increase a template's duration, thereby, com- 
pensating somewhat its lack of the merger and ringdown 
phases. As such, for a given template with a fixed M, 
decreasing its 77 increases its duration. So, relative to a 
target signal's M and 77 values, decreasing a template's 
M or increasing its 77, or both, can produce waveforms 
with durations closer to that of the target signal, in prin- 
ciple. In practice, however, the characteristics of the tar- 
get signals and the modeled templates determine which 
of these three possibilities is realized. 

Note, as well, that for a given target 77, the fitting fac- 
tor recovers after the dip and improves for large values of 
target M since it is easier to fit target signals that have 
a smaller number of cycles. This is consistent with the 
behavior of the fitting-factor as a function of the total- 
mass studied in Ref. HHl. (Similar effects were found 



in Ref. [35j when filtering TaylorEt target signals with 
TaylorTl and TaylorT4 templates.) We do not present 
the results for systematic effects on </>o since it has lit- 
tle astrophysical relevance and is known to incur large 
statistical errors (for SNRs around 10). The error in to 
is more interesting, especially, for triggering searches in 
other electromagnetic observatories, and will be studied 
in Ref. [H. 



C. Analytic approximation using Fisher matrix 



We begin by recalling that Eq. (|2.16[) was obtained 
by dropping quadratic and higher order terms in A9 l 
and, thus, is not expected to be valid for large discrep- 
ancies in the waveform model. The magnitude of errors 
found in our numerical studies makes it manifest that the 
discrepancy between our target signals and templates is 
smaller for smaller M and higher 77. We, therefore, use 
the analytic approximation in Eq. (|2.16l) to compute the 
systematic errors in the parameters in that region, i.e., 
for 7771.2 S [5 - 20]M Q . The results are given for AM/M 
in Fig. [3]and for A77/77 in Fig. [5] (in the right plot). Note 
that we expect the approximation to get worse with in- 
creasing target M and decreasing 77. So the primary pur- 
pose of these figures is to check if the basic trends seen 
in the numerical results discussed above are consistent 
with the theoretical implications of Eq. (J2.16I) . 

Figure [3] nicely corroborates the numerical result pre- 
sented in Fig. [21 namely, that the longer templates (with 
smaller M) provide a better fit for any target signal and 
that the magnitude of the fractional error in the esti- 
mate of M increases with increasing total mass of the 
target signal. For 77, the broad trend of templates with 
larger 77 giving a better fit is visible also in the analyti- 
cal result. However, the specific feature predicted by Eq. 
(I2.16[) that is at variance with the numerical result is the 
behavior of A77/77 near 77 = 0.25. Our numerical results 
show that A77/77 is quite small near 77 = 0.25 and does not 
appear to change even as the target M is increased. Con- 
versely, the analytic formula predicts through the right 
plot in Fig. [5]that the error A77/77 increases with increas- 
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ing M. This apparent disagreement is easily explained 
by the fact that in our numerical simulations we limited 
the templates to have r] < 0.25; no such restriction is 
assumed in the derivation of Eq. (|2.16|) . This is why the 
best fit value of r\ in the right plot in Fig. [5] is larger 
than 0.25 as one approaches the top-right corner, i.e., 
where mi ~ m<x and the target and template waveforms 
start departing from each other. This also suggests that 
allowing for templates with unphysical values of r\ might 
help mitigate some of the loss in SNR arising from in- 
accurate waveform modeling. (See Ref. [HI for possible 
pitfalls of such a strategy.) 

III. DISCUSSION 

In this paper, we studied the loss of SNR and the 
systematic errors in measuring signal parameters ex- 
pected in searches for inspiral-merger-ringdown signals 
from non-spinning BBHs using post-Newtonian tem- 
plates. Our numerical results show that the fitting factor 
can be as low as 60% even for a BBH with m\ = 13M Q 
and m,2 = 20M Q . Also, the estimated total-mass for 
the same pair can be off by as much as 20%. Both of 
these estimates get worse for some higher-mass combi- 
nations. Even the estimation of r/ suffers a nearly 20% 



error for this example, and can be worse than 50% for the 
mass ranges studied here. The implications of the loss of 
SNR and the errors in measuring the masses, such as on 
inferring the effective distance and other astrophysical 
pro perties of the source will be detailed in a future work 
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FIG. 1: The fitting factor obtained from numerical simulations of matched filtering with a bank of 3.5PN TaylorTl templates of 
(complete) inspiral-merger-ringdown waveforms as target signals in AdvLIGO PSD. The target waveforms are parameterized 
by the BBH component masses mi and 7712, each ranging from 13 - 104 Mq. 
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FIG. 2: The fractional error in total-mass (in %) obtained from numerical simulations of the same template bank and target 
signals, in AdvLIGO PSD, as shown in Fig. [1] Target signals with r\ = 0.25 are represented by points along the equal-mass 
line (not shown) extending from the left-bottom corner to the top-right corner of the plot. 




FIG. 3: The fractional error in total-mass (in %) given by the analytic expression Eq. (|2.16[) for AdvLIGO PSD. Above, mi 
and m2 represent the true parameters. 
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FIG. 4: The fractional error in the symmetric mass-ratio (in %) obtained from numerical simulations for the same template 
bank and target signals, in AdvLIGO PSD, as shown in Fig. [T] 




FIG. 5: The fractional error in the symmetric mass- ratio in AdvLIGO PSD (in %). The figure on the left is obtained from 
numerical simulations, and that on the right is computed from the analytic expression in Eq. (I2.16II . 



